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Abstract 

This is the first of a series of papers devoted to the study of classical 
initial-boundary value problems of Dirichlet, Neumann and mixed 
type for the Nonlinear Schrodinger equation on the segment. Consid- 
ering proper periodic discontinuous extensions of the profile, generated 
by suitable point-like sources, we show that the above boundary value 
problems can be rewritten as nonlinear dynamical systems for suit- 
able sets of algebro-geometric spectral data, generalizing the classical 
Dubrovin equations. 

In this paper we consider, as a first illustration of the above method, 
the case of the Dirichlet problem on the segment with zero-boundary 
value at one end, and we show that the corresponding dynamical sys- 
tem for the spectral data can be written as a system of ODEs with 
algebraic right-hand side. 
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1 Introduction 



In the natural phenomena, evolutionary processes very often take place 
in a domain bounded in space and time. At the boundary the system 
has a non-trivial interaction with the external world or experimental 
environment. The mathematical description of such situations is given 
by initial boundary value (IBV) problems for evolution equations. 
Since soliton equations are largely applicable, IBV problems for them 
are of large interest from both mathematical and physical points of 
view. It is the case, for instance, of the IBV problems for the Nonlinear 
Schrodinger (NLS) equation 

iqt = Qxx ± 2g^g, q = q{x, t), q eC (1) 

which describes the complex amplitude modulation of a quasi-monochromatic 
wave packet in strongly dispersive and weakly nonlinear media, and 
then arises in many natural contexts, like in fiber optics, plasma 
physics and fluid dynamics. 

Two IVB problems for soliton equations are well-studied and un- 
derstood. 

1) The initial problem on the line with rapidly decreasing boundary 
conditions at infinity (the Cauchy problem on the line). This prob- 
lem is solved using the scattering transform, first introduced for the 
Korteweg-de Vries (KdV) equation in [21] and later used to solve the 
NLS equation in |46j . In this case the scattering data evolve linearly 
and can be treated as a nonlinear analog of the continuous Fourier 
data. 

2) The IBV problem on the segment [0,L] with periodic boundary 
conditions. Its solution is based on the so-called finite-gap approach, 
first introduced in |i40) . 

More general IBV problems (for example with Dirichlet, Neumann, 
mixed and Robin boundary conditions) for soliton equations are be- 
lieved to be non-integrable. Only for special choices of the boundary 
values the problems remain integrable [T] . |44j and the above two stan- 
dard integration procedures with appropriate reductions on the spec- 
tral data can be applied ( see P for the semi-line case and [Sj for the 
segment one; see also j45j for the semi-line case using the approach 
of IB]); but these cases cover only a limited number of interesting 
physical applications. 

Applying the scattering transform to more general IBV problems 
on the semi- line or on the segment, one meets the following basic 
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difficulty: the time evolution of the scattering matrix depends not 
only on prescribed boundary values, but also on unknown ones. Sev- 
eral strategies have been developed to overcome such difficulty (a list 
of them can be found at the end of this introduction). One of the 
possible strategies, motivated by the classical solutions of the Dirich- 
let and Neumann problems for second order linear PDEs by the sine 
and cosine Fourier transforms, consists in transforming the given IBV 
problem into a certain Cauchy problem on the whole line by introduc- 
ing a proper extension of the profile such that the known boundary 
values are encoded into external point sources and the unknown ones 
disappear from this new formulation. Using this strategy, the Dirich- 
let problem on the semi-line for the NLS equation was investigated in 

More recently, the same strategy was applied to the Dirichlet and 
Neumann problems on the semi-line for the NLS equation in (71 using 
a different formalism. In [Jj the time evolution of the scattering matrix 
^(A, t) of the transformed problem on the whole line is expressed in 
terms of the scattering matrix S{X, t) of the original semi-line problem 
and of the given boundary values. To close the system, one uses the 
non-local connection S{X,t) = S~^{—X,t)aS{X,t) , where o" = 1 for 
the Dirichlet case, and a = for the Neumann case {cr^ is the Pauli 
matrix). Due to the analyticity properties of the scattering matrix 
S{X,t), the above non-local relation can be resolved for S{X,t) using 
the appropriate Riemann-Hilbert problem. 

Motivated by this result, we have recently begun a research activ- 
ity with the purpose of applying the same methodology to the IBV 
problems on the segment [0,L] for the NLS equation. In particular, 
we have found that: 

a) the Dirichlet problem is conveniently reformulated as an odd 
2L-periodic Cauchy problem on the line with 6' sources located at the 
points nL with n G Z. 

b) The Neumann problem is reformulated as an even 2L-periodic 
Cauchy problem on the line with 5 sources located at the points nL. 

c) The mixed problem with, for instance, u{0, t) and Ux{L, t) given, 
is reformulated as an odd 4L-periodic Cauchy problem such that 
q{2L — x,t) = q{x,t), < x < L, with 5' sources located at the 
points 2nL and 5 sources at the points (2n -|- 1)L. Here and later we 
denote by u{x, t) the NLS field in the segment [0, L] and by q{x, t) its 
extension to the whole line. 
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Fig 1. 

Odd 2L-periodic extension for the Dirichlet IBV problem. 



Through these new formulations of the original IBV problems, we 
have been able to obtain the following results. 

1) We have eliminated from the formalism all the unknown bound- 
ary values. 

2) We have transformed the original IBV problem into a Cauchy 
problem for periodic profiles, for which the powerful tools of the finite- 
gap method are at hand. 

3) We have used the algebro-geometric tools of the finite gap 
method to reduce the original IBV problem to a nonlinear system 
of ordinary differential equations describing the time evolution of the 
spectral data of the transformed periodic problem with forcings. This 
system gives a satisfactory analytic description of the IBV problem 
under investigation. 

Since the IBV problem on the segment is equivalent to a periodic 
problem for non-smooth profiles produced by point-like sources, then 
the corresponding spectral theory presents the following novel features 
with respect to the smooth case. 

1) The lengths of the spectral gaps decay slowly in the high-energy 
limit (they are of order 1/ A for the Dirichlet problem and of order 1/ 
for the Neumann one), while, in the smooth case, the gaps lengths 
decay faster than any power of A. 
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2) The number of gaps is always infinite and the finite-gap approx- 
imation converges rather slowly. 

3) The branch points are not conserved any more and their time 
evolution must be included into the set Dubrovin equations. 

4) The evolution of the spectral curve is described in terms of the 
isospectral deformations suggested in 

In this paper we illustrate the above method on the following par- 
ticular Dirichlet problem for the defocussing NLS equation: 

iut = Uxx — ^u^u, u = u{x, t), u G C, < X < L, t > 0. (2) 

n(0, t) = vo{t), u{L,t) = VL{t) = 0, u{x,0) = uo{x), < x < L, 

(3) 

where VQ{t) and uo{x) satisfy the obvious matching conditions v{0) = 
tio(O), uo{L) = 0. The choice of the defocussing case is due to its sim- 
pler analytic structure (the auxiliary spectral problem is self-adjoint) 
and the choice of 0-boundary at x = L simplifies the Dubrovin equa- 
tions. From the point of view of the method used, these simplifications 
are not essential and the solution of the general Dirichlet, Neumann 
and mixed IBV problems, for both the focussing and defocussing NLS, 
will be presented in a more detailed forthcoming paper. 

As it was previously indicated, to study the IBV problem © 
it is convenient to introduce the odd 2L-periodic extension of u{x,t): 

q{x, t)= - 2nL, t)H{x - 2nL)H{{2n + 1)L - x)- 

nez (4) 
u{2nL-x,)H{2nL-x)H{x-{2n-l)L), x G M, t > 0. 

This extension satisfies the following odd 2L-periodic Cauchy problem 
on the line with 5' forcings 

iqt = qxx-2q^q-2vQ{t)^5'{x-2nL), x € M, t > 0, (5) 
q[x, 0) = ^0(2; - 2nL)H{x - 2nL)H{{2n + 1)L - x)- 

nGZ (6) 

uo(2nL - x)H{2nL - x)H{x - (2n - 1)L), 

which will be the subject of investigation of the following sections. We 
show, in particular, that unlike the case of generic forcings, which lead 
to a highly non-local formalism, the point-like sources arising from the 
IBV problem 0, © lead to a system of ODEs for the spectral data 
which can be written in local form after a proper extension of the phase 
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space. This picture is consistent with the mild non-locahty found in 
the semi-hne case [7]. 

The paper is organized as fohows. In §2 we summarize the main re- 
sults concerning the finite gap theory of the defocussing NLS equation 
for smooth periodic profiles. In §3 we present the spectral characteri- 
zation of the discontinuous periodic profiles arising from the Dirichlet 
problem. In §4 we derive the nonlinear system of ordinary differential 
equations describing the time evolution of the spectral data. 

We end this introductory section with a list of other approaches 
to the study of IBV problems for soliton equations, developed during 
the last few years. In [l^I an "elbow scattering" in the {x,t)- plane 
is presented to deal with the semi-line problem for KdV, leading to 
Gel'fand - Levitan - Marchenko formulations. In jl71 118j a different 
approach, based on a simultaneous x-t spectral transform, has been 
introduced and rigorously developed in jl91 l2Uj , to solve IBV problems 
for soliton equations on the semi-line. We were recently informed 
j21j that this method was also extended to the segment case. This 
method allows for a rigorous asymptotics [ini and captures in a natural 
way the known cases of linearizable IBV problems. In (Hj and [Jj 
two alternative approaches to the study of IBV problems for soliton 
equations on the segment and on the semi-line have been presented. 
In the first method, applicable to both the semi-line and segment 
cases, the unknown boundary values are expressed in terms of elements 
of the scattering matrix ^(A,^), thus obtaining a nonlinear integro- 
differential evolution equation for S. In the second method one applies 
a suitable linear operator to StS~^ to eliminate the unknown boundary 
data. It is also shown that the resulting time evolution of S, nonlinear 
and mildly nonlocal, coincides with that obtained applying the point 
source strategy of this paper [7j. Also the formalism presented in [7] 
allows for rigorous asymptotics and captures in a natural way the case 
of linearizable IBV problems. Some integrable boundary conditions 
for soliton equations and their connection with symmetries have been 
investigated in j25j and in references therein quoted. In some non- 
generic cases of soliton equations corresponding to singular dispersion 
relations, like the stimulated Raman scattering (SRS) equations and 
the sine Gordon (SO) equation in light-cone coordinates, the evolution 
equation of the scattering matrix does not contain unknown boundary 
data. The SG equation on the semi-line has been treated using the 
x-t spectral transform fTT^; the SRS and the SG equations on the 
semi-line have also been treated using a more traditional x- transform 
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method respectively in [31^ and in |35j : the x- spectral data used in this 
last approach satisfy a nonlinear evolution equation of Riccati type. 
Apart from the simultaneous x — t transform, all the above approaches 
are based on the traditional 1ST |17j- A different approach, based 
on the Kac-Moody representation for the SG equation in laboratory 
coordinates, was used in ^ to solve the Robin problem on the semi- 
line; it uses critically the finite-speed character of the equation. 

The spectral formalism for studying forced soliton equations has 
been developed by several authors, especially in connection to the 
theory of perturbations (see, for instance, jBHl and ISHl))- Soliton 
equations subjected to one-point-like source were investigated in jT^. 
Soliton equations perturbed by slowly varying forcings were investi- 
gated by averaging procedures in |13j and j31j . We finally remark that 
IBV problems and forced problems for C - integrable equations have 
been considered in [B], 



2 Periodic spectral transform for the 
defocussing NLS 

2.1 Direct periodic spectral transform 

The Lax pair for the defocussing NLS equation: 



iqt = Qxx - 'iq^q, q = q{x, t), 

was found in We use its zero-curvature form: 



(7) 



d^{X,x,t) 



U{X,x,t)'^{X,x,t) 



(8) 



dx 




V{X,x,t)^{X,x,t), 



(9) 



where ^'(A,x,t) is a 2-component vector: 



*(A,x,t) 



^^{X,x,t) 



(10) 



and U{X,x,t), V{X,x,t) are the following 2x2 matrices: 



U{X,x,t) 



iX 
iq{x,t) 



iq{x, t) 
-iX 



(11) 
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iqq qx 
qx -iq<-. 



(12) 



V{\x,t) = 2XU{X,x,t) + 
Equation Q can be rewritten as the following spectral problem 
C^{X,x,t) = X^{X,x,t), (13) 

where 

C=\ (14) 
-q{x,t) id^ \ 

The first symmetry of the NLS hierarchy is the infinitesimal gauge 

qto = iq- (15) 
Its zero-curvature representation has the following form: 



ox Oto 



where 



Vq = -CJ3 



- 

2 ^ 

-- 

^ 2 



(16) 
(17) 



Assume, that our potential is periodic with the period 2L: 

q{x + 2L,t) =q{x,t) (18) 
Denote by T(A, x, a, t) the fundamental solution of 



dT{X, X, a, t) 
dx 

We shall use: 



U{X,x,t)T{X,x,a,t), T{X,a,a,t) 



T{X, b, a, t) = ^(A, 6, t)'^'\X, a, t) 



1 
1 



(19) 



(20) 



where ^{X,x,t) is an arbitrary matrix solution of (jSJ with non-zero 
determinant. 
The matrix 

f{X,t) = f{X,2L,0,t) (21) 

is called monodromy matrix. 

Let us associate spectral data with q{x, t) (we treat JC as an ordi- 
nary differential operator in x depending on the parameter t). 



Lemma 1 The following spectral problems are self-adjoint: 
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1. Main problem: the spectral problem for the operator C on the 
whole line in the space L^(M). 

2. Auxiliary problem: the spectral problem for the operator L on the 
interval [0, 2L] with the boundary conditions: 



The spectrum of the main problem is the union of closed intervals 
of the real line. The complement of this spectrum is the union of 
open intervals called gaps. Let us denote the boundary points of the 
gap number j by and £'2^+1 ) < ^2j+i- The spectrum of the 
auxiliary problem is discrete. Each interval [E2j, E2j+i] (the closure 
of the gap j) contains exactly one point of the auxiliary spectrum. 

For generic potentials the gaps can be enumerated by integers 
— 00 < J < 00, and for sufficiently large j the gap number j is lo- 
cated near the point irj /2L. Each point of the auxiliary spectrum 
is associated with a gap, and we can assume that £?2fc < Afc < E^h+i- 
The lengths of gaps tends to as |j[ ^00. The decay rate depends 
on the smoothness oi q{x, t) and coincides with the decay rate of the 
Fourier coefficients of the potential. 

For non-generic potentials the spectrum of the auxiliary problem 
can also be enumerated so that ■Kk/2L \k\ — >■ oo, but some 

of these points lie inside the spectrum of the main problem. Such 
potentials can be interpreted as the results of "closing" some gaps: 
E2j+i E2j. To construct the inverse transform for the self-adjoint 
problem discussed in our text, only the points of the auxiliary spec- 
trum Aj associated with gaps are essential. 

These spectra possess the following simple characterization in terms 
of the monodromy matrix T(A, t): 

Lemma 2 1. A real number A lies in the spectrum of the main 
problem if and only if \ tr T(A, t)\ < 2. 

2. A number Afe is a point of the auxiliary spectrum if and only if 
the vector 



(0) + -0^ (0) = 0, (2i^) + -0^ C^L) = 0. 



(22) 




is an eigenvector ofT{X,t), or equivalently, 



Tu{X,t)+T2i{X,t) - Ti2{X,t) - T22{X,t) = 0. 
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The Bloch eigenfunction ^'(7, x) is, by definition, the common 
eigenfunction of L and of the shift operator 



/:^(7,j;,t) = A*(7,x,i), ^(7, + 2L, t) = e2*-^^'(^)^'(7, x, t), (23) 



where 7 is a point of a Riemann surface F, which is a two-sheeted 
covering of the A-plane, V denotes the projection from F to the A- 
plane: 'P7 = A and e^'^^'P{i) \^ ^^i. eigenvalue of the monodromy matrix. 
The branch points of F are real and coincide with the ends Ej of the 
gaps. Denote by the oval lying over the interval [£'2^, -E'2j+A:] and 
by a the transposition of sheets. 

The function ^(7) is defined modulo (7r/L)n, n £ Z. It is called 
the quasimomentum. It is odd with respect to the transposition of 
sheets: 

p{a-f) = -p{-f), mod(7r/L). (24) 

The differential of the quasimomentum dp = ^dX is holomorphic 
at the finite part of the spectral curve F. It has exactly 2 zeroes in 
each oval a^: and a^, cra^ 



a 



dp{j, t) 



k ■ 



0, 



lm.p{at) > 0. 



(25) 



(26) 



Denote by Uk the projection of these zeroes to the A-plane. 

Equation (|23() defines ^'(7,2;,*) up to a constant factor. In our 
text we use the following normalization: 



$(7,0,0) = 1, 



where 



^(7, X, t) = 7p\j, X, t) + lP'^{-f, X, t). 

It is also convenient to introduce the following function: 



(27) 
(28) 



E{j,x,t) = {i'\^,x,t) - il;\j,x,t))<!>~\j,x,t). (29) 
By definition: 



^{-f,x,t) 



l+5(7,a:,t) 

1-S(7,a:,t) 
2 



(30) 
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For each x and t, function ^('y,x,t) has exactly one zero 'yk{x,t) 
on the oval Uk- The poles 7^ of $(7,x,i) do not depend on x, t, they 
lie on the ovals over the intervals [E2k, E2j+k] and each oval contains 
exactly one pole 7^, with 7^ = 7^(0,0). Denote by Xk the projections 
of the points 7^ to the A-plane. From the definition it follows that they 
are the points of the auxiliary spectrum corresponding to the gaps. 

Hence the following spectral data are associated with a periodic 
potential of the defocussing NLS equation. 

1. A collection of real branch points -B^ (or, equivalently, the spec- 
tral curve r). 

2. The divisor of poles 7^ in F. 

These data uniquely determine the potential q{x,t). 

2.2 The periodic inverse problem. Finite- gap 
potentials 

A potential q{x, t) is called finite-gap if the surface F is algebraic or, 
equivalently, if the genus g of F is finite. In the finite-gap case we 
define F by the algebraic equation: 

2g+l 

fi^ = RiX), R{X)= lliX-Ek). (31) 

A:=0 

By definition a point 7 of F is a pair of complex numbers 7 = (A, /u) sat- 
isfying (|^T|). Therefore the poles and the zeroes of $(A,x,t) are pairs 
of complex numbers 7^ = {Xk,fJ-k) and complex functions 7fc(x,t) = 
(Afc(x, t), //fc(x, i)) respectively, satisfying pif) . Here < k < g. 

The role of the finite-gap potentials in the soliton theory is analo- 
gous to the role of finite Fourier series in the linear theory. 

In the finite-gap case the function 5(7, x, t) has the following ana- 
lytic properties: 

1. H(7,a;,t) is meromorphic in 7 on the spectral curve F. 

2. H(7, X, t) has at most first-order poles at the divisor points 7fc(x, t) 
and no other singularities. 

3. H(7,x,t) — > ±1 as 7 ^ ±00. 
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Lemma 3 The properties formulated above completely determine the 
function '^{j , X , t) . The explicit formula for '^(j , x , t) is given by: 



E{j,x,t) 



rk{x,t) 



k=0 



(32) 



where 



rk{x,t) 



fik{x,t) 



n (Afc(x,t) -Aj(x,t))' 

j=0,....g 



(33) 



One of the possible procedures for reconstructing potentials by the 
spectral data is the following. The ends Ek of the gaps are constants 
of motion, i.e. they are x, t- independent. The divisor of zeroes 7fc(x, t) 
satisfy the following systems of first-order ordinary differential equa- 
tions in X and t (Dubrovin equations): 



dXk{x,t) 
dx 



29+1 



rk{x,t), (34) 



dXk{x,t) 

at 



-2i 



29+1 



j=0 



A-(x,t) 



rk{x,t). 



(35) 



Solving equations (|34|) . (|35jl for the initial data Afc(0, 0) = A^, //fc(0, 0) = 
Hk, one obtains the divisor of zeroes 7fc(x,t) for all x, t. Using the 
reconstruction formula 



^,*) = E 



k=0 



Xk{x, t) — -{E2k + E2k+l) ^ ^ 



(36) 



one gets the potential q{x,t). 

In the derivation of the Dubrovin equations one assumes that the 
spectral curve has a finite genus. These equations can also be used 
for the infinite-gap case, but it is necessary to be careful with the 
associated convergence problems. 

We use also Dubrovin equations for the infinitesimal gauge (|15() : 



dXk{x,t) 
dto 



-rk{x,t) 



(37) 
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In the finite-gap case, the Dubrovin equations can be exphcitly 
linearized using the Abel transform and the potentials can be written 
in terms of Riemann G-functions associated with the curve T. But we 
do not use this property in our text. 

In the finite-gap case the A-derivative of the quasimomentum dif- 
ferential reads as: 

... n(A-afc) 

^ = '-^ (38) 

OA /i 

Remark. The algebro-geometrical (or finite-gap) solutions of the 
KdV equation were first constructed in jlU], where the zero-curvature 
representation was also first introduced. The complete finite-gap the- 
ory for the periodic Schrodinger operator was developed in ^ ^1 ^2 
UniESlEZl (see the survey article [12] and the books [HI, 0). A pure 
algebraic formulation of the finite-gap integration procedure and its 
generalization for 2+1 dimensional systems were obtained in |3Uj . 

Finite-gap solutions of the defocusing NLS equation were first con- 
structed in 12?]. The first 0- functional formulas for the focussing NLS 
were written in |2H]. A detailed study of the finite-gap NLS solu- 
tions can be found in An interpretation of generic periodic 1- 
d Schrodinger potentials in terms of 0-functions of infinitely many 
variables was suggested in f3H^. A generalization of this approach to 
generic matrix ordinary differential operators was suggested in |43j . 
The spectral transform for the operator (jSJ with potentials in the 
Hilbert space L^([0,2L]) was developed in PU] . 

In the finite-gap case we mark a point of 7 € F by a pair of complex 
numbers 7 = (A, ^). But, to study the limit g — > 00, it is convenient 
to renormalize /x. Assume that we have a point inside the k-th gap: 
A G M, < A < E2k+i- Then we define a new variable fl by 



n (^--^j) 



j#2fe 



sgn|n(A-A,)| , (39) 



A' = (^2fc+i-A)(A-S2fe). (40) 
(// is pure imaginary inside gaps and we choose fl to be real, therefore 
we have the constant i in (|39|) ). 

Since the residues are pure imaginary, it is convenient to intro- 
duce the real quantities ffc(x,t) as follows: 

rkix,t) =2ifkix,t). (41) 
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From it follows: 



\J (Afc — E2j){Xk — E2j+i) 
l(Afc-A,)| 



(42) 



2.3 Isoperiodic deformations of the spectral 
curve 

Finite-gap solutions corresponding to generic surfaces are quasiperi- 
odic in x and t. To construct x-periodic solutions with a prescribed 
period, it is necessary to put nontrivial additional constraints on the 
Riemann surface. These constraints for the KdV equation were formu- 
lated in |36j in terms of conformal maps. After minor modifications, 
the approach of [HHI can be applied to the defocussing NLS equation 
but, for our purposes, an alternative approach, suggested in [21], is 
more convenient. The main idea of is to use the so-called isoperi- 
odic deformations, i.e. ordinary differential equations on the moduli 
space of Riemann surfaces such that all x quasiperiods of the potential 
are conservation laws. Therefore, if we have one periodic potential, 
we can construct new solutions integrating these deformations. The 
starting point can be chosen in the neighborhood of the zero potential 
using the perturbation theory. 

Isoperiodic deformations arose in [,V2\ . [4^\ in the form of 

ODE's with a rather complicated transcendental right-hand side. In 
|24j it was shown that, by extending the phase space, these flows can 
be transformed to a very simple form with a rational right-hand side. 

Consider the following collection of g + 1 flows of the branch points 
Ej and of the projections aj of the quasimomentum differential zeroes 
E2j < OLj < E2j+i: 



dEj Cfc 



with j = 0, ...,2s- + 1, k = 0,...,g; 



(43) 



dTk Ej - Ofc 



daj _ Cfc 



with j = 0, 



k-l,k + l,...g; 



(44) 



dTk oij - afc 




(45) 
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where 



2g+l 
n {(^k-E^] 
j=0 



\[{ak-ajY 



(46) 



Here denotes the parameter of the flow associated with the point 

ak- 
in the fohowing properties of equations (|1H )) -H45 |) were proved: 



Lemma 4 1. These flows preserve all x-quasiperiods of the corre- 
sponding NLS equation solutions. In particular, if the starting 
solution is x-periodic with the period 2L, it remains 2L-periodic 
in X. 

2. These flows commute pairwise. 

3. Assume that the parameters are normalized by the following 
condition: they tend to if the lengths of all gaps tend to 0. 
Then 

) = =Fi7-fe, mod{-K/L). (47) 
4-. The variation of the quasimomentum by these flows is given by 
5p(7, f) Ck 5p(7, r) 



dTk 



A 



dX 



(48) 



3 Asymptotics of the spectrum for dis- 
continuous periodic profiles 

In the Fourier theory for periodic smooth functions it is weh-known 
that the coefficients of the series decay faster than any power of k. 
If the function /(x) is piecewise smooth and discontinuous at a finite 
number of points xi, . . . , x„ , then the coefficients decay as 1/A; and 
the leading terms are: 



j=i 



(49) 



Here 2L is the period of the function and 5j is the jump at the point 
^j- ~ fi^j') ~ f(^j)^ where f{x^) = lim /(x ± |e|), e — > 0. 
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In addition, if the right and left derivatives of /(x) coincide at ah 
discontinuity points, then the corrections are of order 

The periodic spectral transform can be treated as a nonlinear ana- 
log of the discrete Fourier transform. In the small potential limit the 
lengths of the gaps are |2/fc| and the phase of fk is encoded in the 
position of the divisor. 

Since the Dirichlet IBV problem is reformulated as an odd 2L- 
periodic Cauchy problem with jumps 2vo{t) at the points 2nL and 
jumps —2vL{t) at the points (2n + l)L, then we have to find a spectral 
characterization of this kind of discontinuous profiles. 

Proposition 1 Let the potential q{x) be odd, 2L-periodic and suffi- 
ciently smooth outside the discontinuity points uL and have the jumps 
2vQ{t), —2vL{t) at the points 2nL, {2n + 1)L respectively. Then the 
odd character of q{x) implies the following symmetry of the spectrum: 

Ek = -E^k, Afc = -A_fc, keZ (50) 

and the above discontinuities imply the following asymptotics: 




where 

dk = dk{t) = vo{t) -{-!)'' V Lit), 

hit) = J dx\uix,t)\^ ^^^^ 



and voit) = uiO,t), v^it) = uiL,t). 

The absence of 0(1/A;^) corrections, due to the symmetry (|5U|) . will 
be essential for us. 

The idea of the proof is the following: in each continuity interval 
equation Q can be asymptotically diagonalized by the gauge trans- 
formation 

^(A,x) = G(A,x)^i(A,3;), (53) 

where G(A, x) is an asymptotic series in 1/A. All coefficients of G(A, x) 
can be computed explicitly in terms of qix) and its derivatives. Taking 
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into account the discontinuities at 0, L and 2L, one gets the fohowing 
representation for the transition matrix in the interval [0, 2L]: 

T(A) = G(A,0-)e^^P2(^)"^G-i(A,L+)G(A,L-)e*^P^(^)"^G-i(A,0+), 

(54) 

where pi,2(A) are the averages on the intervals [0,L] and [L,2L] re- 
spectively of the first component of the diagonalized U matrix. 
By Lemma the branch points are determined by the equation 

trr(A) = ±2 

and the divisor at x = is determined by the equation 

Til (A) + T2i(A) - ri2(A) - T22(A) = 0. 

The asymptotics (|?T|) can be obtained expanding these equations in 
inverse powers of A. 

Lemma 5 Let T be an infinite genus spectral curve with real branch 
points Ek such, that 



ttA; ^ / 1 
E2k = — + O - 



tt/c ^ / 1 



(55) 



and each oval Ofc contains exactly one divisor point 7^. Then the 
quantities r^ and r^ given by and o,re well-defined. 

The proof is rather straightforward. 







Re q(x) 




^"^^^^^ Im ( 


^(x) 







Fig 2. 

This potential has 80 gaps (it corresponds to 40 Fourier harmonics). 
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4 Time evolution of the spectral data 



In this section we derive the equations describing the time evolution 
of the spectral data - the spectral curve and the divisor - for the 
Dirichlet problem with zero at one end 



u{0,t) = vo{t), u{L,t) = VL{t) = 0, u{x,0) = uo{x), < x < L. 

(56) 

Let 

vo{t) = A{t)e'^('\ (57) 

where A{t) and B{t) are real functions. 

The calculations show that the evolution of the divisor is described 
by the sum of two separately divergent vector fields. But if the bound- 
ary function vo{t) is real, both divergences simultaneously vanish. This 
suggests to use the following gauge: 

g(x,t) =^(x,t)e*^W (58) 

and to calculate the Dubrovin equations for the function q{x,t). It 
satisfies the following equation: 

iqt{x, t) = qxx{x, t)-2\q{x, t)\'^q{x, t)+B{t)q{x, t)-2A{t) ^ 5'{x-2Ln). 

(59) 

The right-hand side of (|59|) is a sum of 3 vector fields: 

^=Dsq + Dgq + Dfq, (60) 

where 

Dsq = -iq{x,t)xx + '^i\q{x,t)\'^q{x,t), (61) 
Dgq = -iB{t)q{x,t), (62) 

Dfq = 2iA{t)'^6'ix-2Ln). (63) 

Starting from now we denote by Xk{x,t) the divisor corresponding 
to the gauged potential q{x,t). It is natural to write the evolution 
of the branch points, of the zeroes of the quasimomentum differential 
and of the divisor as the sum of these 3 flows: 

= DsEkit) + DgEkit) + DfEkit), (64) 
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dt 

dt 



= Dsak{t)+Dgak{t) + Dfak{t), (65) 
Ds\k{^,t) + Dg\k{<d,t) + Df\k{^,t). (66) 



The action of the vector fields Dg and Dg can be obtained as the 
g ^ oo hmit of the finite-gap formulas, but it is necessary to check 
the convergence of these expression. From the results mentioned in 
Section 2 and, in particular, from equations (jSSl), ((SZ|)) we have 



DsEkit) = DgEk{t) = D,ak{t) = DgUkit) = 0, 



DsXk{0,t)=4 



DgXkiO,t) 



dBit) , , 



rfc(0,t), 



(67) 
(68) 

(69) 



Lemma 6 Let the spectral data have the asymptotics i51\) of Propo- 
sition^ with dk{t) = A(t) G M. Then the right-hand side of i68\) is 
well-defined. 

The proof is straightforward. 

The rest of the section is devoted to the proof of the following Main 
Theorem, describing the effect of the forcing on the dynamics of the 
spectral data. 

Theorem 1 The effect of the forcing on the spectral data is described 
by the following vector fields: 



Cfc(t) 



DMt), 



+ 



otjit) - ak{t) 

2^ Ek{t)- a J (t) ^ ak (t) - aj (t) 



Cj{t) 



DfXkiO,t)= 4A(t)Afc(0,t)^^^||±^+2A(t)2 
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DfTkit) + 
DfT,{t), 

fk{0,t), 



(70) 
(71) 



(72) 



where 



DfTkit) 



-2afc(t)^^^^l-H(a+,0,i)H(a^,0,t) 



L 



•{a+,0,t)-E{a;,0,t) 



^(7,x,t) 



A - \k{x,t) 



J--'- \X-Xj{x,t)\ ^ 



(73) 
(74) 
fk{x,t) 



{X-Xk{x,t)) 



(75) 

and Cfc are defined by jj^A )- 

These vector fields can be transformed into an algebraic form by 
introducing the following extended phase space: {Ek, ak, Xk, f^k, M^} ■ 
The equations describing the time evolution of the remaining parame- 
ters have the following form: 



^/^.(o,t) = ^^ 



Xk{Q,t) - jrE2k{t) 4fXk{0,t) 



dt 



Xk{0,t) - E2k{t) 



+ 



di^2fc+l(t) 



-Mkit) = iL{Mk{t)+l)ll 



Xk{0,t) - E2k+i{t) 

(76)" 

\{Xk{0,t)-a,m 
j^k VMO,t) - E2jmXk{0,t) - E2j+i{t)) 

(77) 



XkiO,t) - akjt) 
fj'kiO,t) 



|Afc(0,t) 



E 



f Xk{0,t) - a^it) 



Dnit) 



Remark. The above dynamical system contains infinite sums and 
products. Their convergence follows from direct calculations using the 
asymptotic formulas (|51j) . 

The proof of Theorem ^ is based on the following formula: 



2L 



5T{X,t) = J T{X,2L,x,t)6U{X,x,t)T{X,x,0,t)dx. (78) 



Let us introduce the following notations: 

r(7, t, to) = §^^7, 0, to)f{X, t)§(7, 0, to), (79) 

where ^(7, x, t) is the following 2x2 matrix solution of Q combined 
from the Bloch solutions: 

^{^,x,t)=Mj,x,tMaj,x,t)] = ( ^;(-^'-';) 

V V' (7,a;,t) 'il^'{aj,x,t) 

(80) 
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By definition, 



„2iLp(7,to) Q 

nX,to,to)={ Q g-2.Lp(,,o) )• (81) 



A direct calculation shows that 

2L 



6fiX,t,to) =f{X,to,to) [ ^~\X,x,to)5UiX,x,to)^i\,x,to)dx. 

t=to J 


(82) 

For the variation of the potential q defined by (|63() we obtain 

2L 

t=to 'W J V -V' (7,a;,to) V (7>2;,to) / 



i5\ -2A{tQ)[5'{5:) + 5'{x-2L)]5f 



2A{tQ)[6'{x) + 5'{x-2L)]5f -i5\ 

(83) 

^^^(7,5,^0) ^'^i(^7,x,to) 
ip'^{'y,x,to) V^((T7,x,to) 

_ 1 f e^'^P(^''o) \ / An A12 V ^ . 

~W\ e-2*Wo) J [ A21 A22 

where W is the Wronskian determinant 

W = 1/;^ (7, X, to)V'^(f^7, X, to) - ^^^(7, x, to)iP^ {a-f, x, to), (84) 

2L 

An = i5X [ '>p'^{'y,x,to)Tp'^{crj,x,to) +ij'^{'y,x,to)ip^i(7'y,x,to)dx + 

(85) 







2L 



+2^(to)5/ / (V'^(7,x,io)V'^(f^7,2;,to) -V'^(7,a;,to)^^(fT7,2;,to)) [S'ix)+6'ix-2L)]dx, 



2L 



A21 = -2i5A / V^(7,2;,to)V'^(7,a;,to)da;+ (86) 



2L 



+2A{to)5f J ((V'(7,x,to))'-(^'(7,^,io))') [6' {x) + 6' (x - 2L)]da 
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(we do not use the formulas for A12, A22 in our calculations). 

Let us first calculate the deformation of the spectral curve. Prom 
it follows that 



dt 



dp{a+{t)),t) 



t=to 



dt 



dpi'j, t) 



t=to 



dt 



t = to 



Then, from ()83() . it follows that 



dn 

dt 



A 
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t=to 



2LW 



with 5\ = 0. 

The quadratic combination of eigenfunctions in ()85() is periodic, 
therefore we have 

An = 2^(to) X (87) 



X / ['ip^{al,x,tQ)^^{a^,x,tQ) -i)^{al,x,tQ)^'^{a^,x,tQ))5'{x)dx 



= -2A{tQ) U^{al,x,tQ)'il?-{a^ ,x,to) - ip'^ia'l: ,x,to)^'^{af^ ,x,to)]' 

= -Aiakito)A{to) [iPHa+,0, to)i^\a^,0, to) + V'^(a+, 0, to)ij\a^,0, to)] 

Using the fact that the Wronskian does not depend on x, the defi- 
nition (|5n|) and the arbitrariness of to, we obtain equation (|74() . Using 
the isoperiodic deformation formulas (|43j) - (|44)) . we finally obtain the 
evolution equations ((70)) . ((7T|) of the spectral curve. 

Let us calculate now the evolution of the divisor with respect to 
the vector field D j . Combining Lemma |21 with the definition of T 
we obtain the following property of the divisor points A^: the vector 

remains the eigenvector of T after variation, or equivalently. 







A2i=0. Therefore 



where 



DfXk{0,t)=iA{t)^, 



2L 



AT 



{ii;\jkiO, t),x,t)f - (V^^(7fe(0,t),x,t))2) [6'ix)+6'{x-2L)]dx 
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2L 

J {{i;\jk{0,t),x,t)f - {i;\^k{0,t),x,t)f)' [5{x) + 5{x - 2L)]dx, 



2L 

2^= / i^H7kiO,t),x,t)i;\jkiO,t),x,t)dx. 



We used here the definiton of the divisor points 7^ by the condition 
ip^{'yk{x,t),x,t) + ip'^{'~fk{x,t),x,t) = 0. It is convenient to normahze 
the function ^'(7,a;,t) near the point 7^(0, t) by assuming 

V^i(7,0,t)-V'(7,0,t) = 2. (89) 

As a corollary we obtain 

^i(7fc(0,t),0,t) = 1, V'(7fe(0,i),0,t) =-1, (90) 

V^i(7,0,t)+V'(7,0,t) = — ^TTY- (91) 
Let us calculate M: 

{{^\^,x,t)f-{^\^,x,t)f)' = 

= -2i\{',jj^{^, X, t))2 + (^i(^, X, t) f)-2i{q{x, t)+q{x,t))i;\^, x, t)tp\-f, x, t); 
therefore 

2L 

M=-2iXk{0,t) I ((^2(^^(o,t),x,t))2 + (V'i(7fc(0,t),x,t))2) [5{x)+5{x-2L)]dx- 



(92) 

2L 

-2i / (^(x, t)+q{x, t))V'^(7fc(0, t), X, t)V'^(7fc(0, t), x, t)[6{x)+5{x-2L)]dx. 







If /(x) is a continuous function, then 



e 



j f{x)5{x)dx = ^f{0), (93) 


where e is any positive constant. Therefore the first integral in (|92|) is 
well-defined: 

2L 

((^2(^^(0, t),x, t)f + (Vn7fc(0, t), X, t)f) [6{x) + 5(x - 2L)]dx = 

(94) 
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,4iLp(7fc(0,t)) 



+ 1. 



But the second integral depends on the regularization of the (5-function. 
In our problem the (5-term arises as the compensation of the disconti- 
nuity in q{x, t). Therefore it is natural to write, near the point x = 0, 



6{x) 



2A{t) 



Re qx{x, t) + regular part. 



We use also the fact that the function Iuiq(x,t) is continuous and 
vanishes at the point x = 0; therefore 



2L 



2(Re q{x, t))^^ (7^(0, t),x, t)V'(7fc(0, t),x, t) [5{x) + 5{x - 2L)] 



lim 



2L 



+ 

.0 2L 



A{t) 



Re q{x, t) Re qx{x, t)V'^(7fc(0, t),x, t)ip^ {'jk{0, t), x, 



(95) 



Combining and (|95() we obtain: 



Af= -2iAfc(0,i) 



j4iLp(7fc{0,i)) , ^ 



iA{t)[e 



4iLp(7fc(0,t)) 



1]. (96) 



Let us now calculate P: 

2L 

T> = lim — 

7'-+7fc(0,t) 2 



— lim 

2 7'^7fc{0,t) 



d 



^^^(7^(0, t),x, t)ij\-f', X, t)+ip\-/', x, t)ilj\-fk{0, t), X, t)dx 
^'^7^(0, t),x, t)i,^ (7', X, t) - V'HV , X, t)V''(7fc(0, t),x, t) 



A'-Afc(0,t) 

„4iLp(7fc{0,i)) _ 



2L 




X ^ [V'' (7fc (0, t) , 0, t)V'' (7', 0, t) - ^1 (7', 0, t)i;^ (7fe (0, t),0, t)] 



7'=7fc(0,t) 



,4iLp(7fc(0,i)) _ ^ 



d 



— [^\J',0,t)+^^J\i,0,t)] 



7'=7fc{0,i) 
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Taking into account we obtain 



V 



„4iLp(7fe(0,t)) 



r-fc(0,t) 



(97) 



and, finally, equation ((7^ . 

Equation (|76|) can be obtained by direct differentiation of ()40|) and 
formula ()77|) follows immediately from H48() and (|73() . 



5 Concluding remarks 

1. In our method we study the IBV problem 0, © through the 
following steps: 

(a) From the initial condition uq{x) we build its odd periodic 
extension g(a;, 0) as in © and then we construct, through 
the direct periodic spectral transform summarized in Sec- 
tion the initial spectral data {£^^(0), afc(O), Afc(0, 0), 
^fc(0,0),Mfc(0)}. 

(b) Solving the system of ODEs (|S11) - (d, dZHl), dZZl) with 
the above initial conditions, we obtain the evoluted spectral 
data {Efc(t), afc(t), Afc(0, t), /Xfc(0, t),Mk{t)}. 

(c) Solving the x-Dubrovin equations pijl we obtain, from (0, t) , 
Afe(x,t). 

(d) Using formula (fS^ we reconstruct the potential q{x^ t), whose 
restriction to the segment [0, L] gives the solution u{x, t) of 
the IBV problem @, ©. 

(e) It is well-known that the x-Dubrovin equations can be lin- 
earized, at fixed t, on the Jacobi variety constructed from 
the spectral data {Ek{t),^kiO,t)}. Correspondingly, the so- 
lution q{x,t) can be written explicitly in terms of Riemann 
^-functions. This is an alternative procedure to the above 
steps ((Tc)l , (jldp . However we remark that this Jacobi variety 
depends on t in a highly non-trivial way, in contrast to the 
integrable cases. 

2. Althouh the x-Dubrovin equations can be linearized on the Ja- 
cobi variety and the solutions can be written down explictly in 
terms of Riemann ^-functions, it turns out that, from the nu- 
merical point of view, it is often (and, in particular, in our case) 
more convenient to integrate the Dubrovin equations directly. 
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This idea was communicated to one of the authors (P.G.) by 
A.Osborne long times ago, and it was used by the author in his 
numerical experiments. 

3. We believe that our approach be the proper nonlinear analog of 
the discrete sine-Fourier transform method, which allows one 
to solve the Dirichlet problem on the segment for the linear 
Schrodinger equation. 

4. We consider this approach - in which the original IBV problem 
is reformulated as a forced periodic IBV problem that, with- 
out forcing, is known to be integrable - as a convenient way 
to treat the original IBV problem. One of the reasons is the 
following. It is well-known that, to analyse a perturbed (by a 
forcing) nonlinear equation, one has to calculate the spectral de- 
composition of the forcing in terms of the eigenfunctions of the 
linearized equation. But, for soliton equations, this decomposi- 
tion can be explicitly written in terms of the auxiliary spectral 
problem. Therefore the calculation of the time evolution of the 
spectral data due to the forcing can be considered as the proper 
nonlinear analog of the Fourier analysis. 

5. A detailed study of the system of ODEs - (f7^ . (fTH)) . (f77|) 
will be the subject of future investigation. Preliminary numer- 
ical experiments seem to show that, starting with zero initial 
condition and applying a smooth localized forcing, only a finite 
number of nonlinear modes remain essentially non-zero when the 
external forcing vanishes. Correspondingly, only a finite subset 
of the above system of ODEs is essential for gaining a correct 
picture of the process. 

6. The method presented in this paper can obviously be applied to 
study the general Dirichlet, Neumann and mixed problems for 
the focussing and defocussing NLS equations on the segment. 
These general results will be presented in subsequent works. 
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